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If the inflationary era is preceded by a radiation dominated era in which the inflaton too was 
in thermal equilibrium at some very early time then the CMB data places an upper bound on 
the comoving temperature of the (decoupled) inflaton quanta. In addition, if one considers models 
of “just enough” inflation, where the number of e-foldings of inflation is just enough to solve the 
horizon and flatness problems, then we get a lower bound on the Hubble parameter during inflation, 
Hi n f, which is in severe conflict with the upper bound from tensor perturbations. Alternatively, 
imposing the upper bound on Hi n f implies that such scenarios are compatible with the data only if 
the number of relativistic degrees of freedom in the thermal bath in the pre-inflationary Universe 
is extremely large (greater than 10 9 or 10 11 ). We are not aware of scenarios in which this can be 
satisfied. 


I. INTRODUCTION 

The era in the early history of the Universe before cos¬ 
mic inflation is completely unknown. It has been ar¬ 
gued [1-3] that before the era of observable inflation, 
the Universe could have been sitting at a local minimum 
(metastable vacuum) in the effective potential of the un¬ 
derlying microscopic theory and then have tunnelled (via 
bubble nulceation) to a sufficiently flat direction causing 
slow-roll inflation. The standard thermal history of the 
Universe begins after the end of inflation when the pro¬ 
cess usually referred to as reheating causes the Universe 
to get into a thermal state. 

In this work, we consider a different picture of the early 
Universe. In our picture, the pre-inflationary early Uni¬ 
verse is in a thermal radiation dominated state but as the 
Universe expands, the density of radiation decreases and 
eventually falls below the energy density corresponding 
to the rolling field just like the energy density of the Uni¬ 
verse gets dominated by a non-zero cosmological constant 
or quintessence field at late times. We ask what can the 
present cosmic microwave background (CMB) data tell 
us about this possibility. 

One crucial consequence of a pre-inflationary thermal 
history is that the quanta of the inflaton fluctuations 
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S(f> (the excitations about the classical inflaton back¬ 
ground) could have been in thermal equilibrium very 
early in the pre-inflationary radiation era. If so, it 
would then subsequently decouple but would retain a 
thermal Bose-Einstein distribution even after decoupling. 
Hence the held 5<f> during inflation would be in a ther¬ 
mal state ((a^ak'} = [exp(fc/T) — l] _1 ^ 3 (k — k'), where 
T is the inflaton comoving temperature), rather than in 
the Bunch-Davies vacuum state. This effect was studied 
in Ref. [4] where it was found that the primordial scalar 
power spectrum gets an additional factor of coth[fc/(2T)] 
and the CMB temperature anisotropies get enhanced at 
large angular scales by an amount which depends upon 
the comoving temperature of the inflaton quanta dur¬ 
ing inflation 1 . From the measurements of temperature 
anistropies of the CMB, one can put an upper limit on 
the comoving temperature of the inflaton. This then has 
implications for the duration of inflation [4]. A more 
detailed analysis, using COSMOMC, and including the 
modification of the inflaton mode functions due to a pre¬ 
inhationary radiation era, has been carried out by us in 
Ref. [8]. However in Ref. [8] we treat both the comoving 
temperature T and the number of e-foldings of inflation 
more than that required to solve the horizon and flatness 
problems (SN) as independent parameters. We obtain 


1 In Refs. [5], [6] and [7] a thermal state of decoupled gravitons 
was considered which enhances the temperature anisotropy and 
the B mode polarisation angular power spectrum of the CMB at 
large angles. 
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a stronger constraint on the comoving temperature than 
in Ref. [4], Furthermore we find that the best fit value 
of the number of e-foldings of inflation is very close to 
the minimum duration of inflation required to solve the 
horizon and flatness problem (though the improvement 
in likelihood is not very much over the standard power 
law spectrum). 

Other previous studies of a scenario with a pre¬ 
inflationary radiation dominated era have also shown 
that the CMB data favours scenarios with just enough 
inflation to solve the horizon and flatness problems (be¬ 
cause of the low power on large scales in the TT angular 
power spectrum) [9-16] (see however Ref. [17]). We point 
out that there also exist good theoretical reasons for con¬ 
sidering models of “just enough inflation” in which SN is 
small. By considering anthropic bounds on the value of 
the curvature and including certain statistical arguments 
on the parameters of inflation, it is argued in Ref. [1] 
that the number of e-foldings of inflation should not be 
much larger than the minimum required by observations. 
In Ref. [18] it is shown that by considering a certain 
natural canonical measure on the space of all classical 
universes the probability for N e-foldings of inflation is 
suppressed by a factor of exp(—31V). Moreover, more in¬ 
flation means that the inflaton potential is required to 
stay flat for a huge range of field values and can require 
super-Planckian excursions. Therefore below we shall a 
consider a scenario of just enough inflation. 

None of the studies in Refs. [9-17] consider the pos¬ 
sibility that the inflaton too could have been in thermal 
equilibrium during the pre-inflationary radiation era, as 
has been discussed above. In the scenario of Ref. [8] 
we presume that the inflaton was in thermal equilibrium 
at some early time in the pre-inflationary radiation era 
and that it subsequently went out of equilibrium and de¬ 
coupled, and so include a (frozen) thermal Bose-Einstein 
distribution for the quanta of the inflaton field. We be¬ 
lieve this is a natural extension of scenarios with a pre- 
inflationary radiation dominated era. After all, the infla¬ 
ton must have some couplings to other fields - some are 
needed for sufficient reheating to occur after inflation. 

Our current analysis below shows that such a scenario 
with just enough inflation and a frozen thermal distribu¬ 
tion for the inflaton quanta is compatible with the CMB 
data only if the effective number of relativistic degrees of 
freedom for the pre-inflationary radiation plasma is ex¬ 
tremely large, in fact, larger than 10 9 or 10 11 , depending 
on whether or not any entropy is transferred to the ra¬ 
diation after the inflaton decouples. This had not been 
recognised in Ref. [8]. 


II. THE MODEL 

Previous studies of pre-inflationary radiation era had 
either considered modifications of the mode functions of 
the inflationary scalar perturbations due to presence of 
the prior radiation era [9-17, 19, 20], or a thermal distri¬ 


bution of inflaton fluctuations associated with the infla¬ 
ton being in thermal equilibrium at some very early phase 
in the radiation era [ ]. While the first contribution low¬ 
ers the quadrupole moment of the CMB TT anisotropy 
spectrum, the other contribution enhances the power at 
large angular scales. We, in Ref. [8], included both these 
effects as a consistent approach to determining the sig¬ 
natures of a pre-inflationary radiation era. We obtained 
the scalar power spectrum in Ref. [8] as 

V n (k) = A (J^j coth C(k), (1) 

where C(k) reflects the modification of the inflaton mode 
functions and hence the power spectrum due to a pre- 
inflationary radiation era. The modified mode functions 
are obtained by setting Bunch-Davies initial conditions 
on the inflaton field in the pre-inflationary radiation dom¬ 
inated era, and matching solutions for the field fluctua¬ 
tions and their time derivatives at the (instantaneous) 
transition from the radiation to the inflationary era (as 
opposed to simply setting the Bunch-Davies initial condi¬ 
tions in the inflationary era). Denoting the physical tem¬ 
perature of the inflaton by T( t ), the comoving tempera¬ 
ture T = a(t)T(t), with a = 1 today. The coth (^) fac¬ 
tor appears due to the thermal distribution of the inflaton 
fluctuations: l+2rifc = l+2[exp(fc/T) — l]” 1 = coth (^). 
If one considers only the modification of the inflationary 
mode function due to the presence of a pre-inflationary 
radiation era [10-17, 19, 20] then the power spectrum 
would not contain the coth[/c/(2T)] factor, and similarly 
considering only the thermal distribution of the inflaton 
quanta the power spectrum would not contain the C(k ) 
factor [4]. 

We carried out a detailed COSMOMC analysis in Ref. 
[8], using Planck 2013 and WMAP nine year data and 
the power spectrum in Eq. (1). In addition to the 6 
standard parameters for the Monte Carlo simulation we 
included two extra parameters - SN , the number of e- 
foldings of inflation in excess of the minimum required to 
solve the horizon and flatness problems, and T, the co¬ 
moving temperature. (The factor C(k ) in the power spec¬ 
trum depends on SN through the time when the change 
from the radiation era to the inflationary era takes place.) 
For an inflationary scale of 10 15 GeV we found that the 
best fit value of SN is 0.08 (the marginalized value is 
consistent with zero with 68 % confidence). This indi¬ 
cates that the data favours a ‘just enough’ inflationary 
scenario. As discussed earlier, there are also theoretical 
motivations for considering such a scenario. We also ob¬ 
tained an upper limit on the comoving temperature as 
T < 1.3 x 10 -4 Mpc -1 at 68% C.L. (for the GUT and 
the electroweak scale). 

We will now show that this scenario with a frozen ther¬ 
mal distribution of inflaton fluctuations and ‘just enough’ 
inflation preceded by a radiation dominated era faces a 
severe conflict as it either implies a very large Hubble 
parameter during inflation that conflicts with the known 
upper bound, or requires an extremely large number of 
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relativistic degrees of freedom (greater than 10 9 or 10 11 ) 
in the early Universe. 


III. ANALYSIS 


In the following analysis we use the subscript 0 to re¬ 
fer to the present epoch. Let t* be the epoch when the 
mode k 0 corresponding to our current Hubble radius was 
crossing the Hubble radius during inflation. Then 

fco = a*H * = a 0 H 0 , (2) 

where a* = a(f*), H * = if(f*), H(t 0 ) = H 0 and 
do = 1. Presuming that the Hubble parameter is approx¬ 
imately constant during inflation, we set ff* = H ln f = 
M? ni /(y/3Mp), where M- m i and Mp = 2.4 x 10 18 GeV are 
the energy scale of inflation and the reduced Planck mass 
respectively. Then 


% = T/a » 
TH ini 
H 0 


(3) 


If there were SN e-foldings of inflation before t* then 


Ti = % exp SN , 


(4) 


where 7 1 is the physical inflaton temperature at the be¬ 
ginning of inflation. 

We assume that the radiation energy density at the 
onset of inflation at U is equal to that of the inflaton, i.e. 
Pr = P<p- Then the radiation temperature at U 

T 7 = (-2— (5) 

\7r 2 ffrel/ 

/ 90 \ -- 

= -- \JH in (Mp . ( 6 ) 

V 71 " Wreiy 


As mentioned above for GUT scale inflation the best fit 
value of SN is 0.08. Assuming the energy scale to be 
much smaller (e.g. the electroweak scale) implies that 
the best fit SN is 0.02. If we assume that SN « 0.05 and 
T < 10 _4 Mpc _1 , and take H 0 ss (4400Mpc) -1 [22] and 
g re i ~ 100, then Eq. (8) implies that 


or, 


H inl 

Mp 


> 1 , 


( 10 ) 


Mint > Mp . (11) 

But from the upper bound on r at k = 0.002 Mpc -1 
(v.i.z. r < 0.1), we have [23] 


Hint 

Mp 


< 3.7 x 10 -5 , 


( 12 ) 


or Mi n f < 1.9 x 10 16 GeV. Clearly the lower bound on 
H[nt in Eq. (10) is in severe conflict with the CMB con¬ 
straint. (Note that for k = 0.002 Mpc -1 the coth factor 
and C in Eq. (1) reduce to ~ 1 and hence we can use Eq. 
(12)0 

As is clear from Eq. (8) this result is a consequence of 
picking a small value of SN (larger values of SN > 6 will 
decrease the lower bound on Hint in Eq. (10) to accept¬ 
able levels). Our result is thus a serious constraint on 
inflationary models with just enough inflation which are 
preceded with a pre-inflationary radiation era where one 
assumes that the inflaton too could have been in thermal 
equilibrium at some early time prior to inflation. 

We now study possible solutions to this conundrum. 
Eq. (8) implies that Eq. (10) can be recast as 


Hint 10 

Mp y/g ie \ 


(13) 


Here g re i is the effective number of relativistic degrees of 
freedom for the plasma filling the Universe at the on¬ 
set of inflation. If the underlying microscopic theory 
is the Standard Model of elementary particle physics, 
g re i = 106.75 while if it is MSSM, then g rel = 228.75 [21]. 
If the number of relativistic degrees of freedom stays the 
same between the epoch of decoupling of the inflaton and 
the epoch of onset of inflation (i.e. no species annihilates 
or becomes non-relativistic in between), we can assume 
that the inflaton temperature equals the radiation tem¬ 
perature at the onset of inflation, i.e. 


Ti=Ty- 


(7) 


Then combining Eqs. (3), (4), (6) and (7), we get 


Hint _ ( 90 V /2 (Ho\ 2 2SN 
Mp \n' 2 g re ]J \T ) 

Below we shall take 

T < lx 10 _4 Mpc” 1 . 


( 8 ) 

(9) 


However, Eq. (12) then implies that 

ffrel > 10 11 . (14) 

This is an extremely large number of relativistic degrees 
of freedom in the early Universe. 

At this point it is worthwhile recalling another as¬ 
sumption we have made in arriving at our result. As 
Eq. (7) indicates, we have assumed that the temper¬ 
ature of the decoupled inflaton quanta and that of the 
pre-inflationary radiation plasma is the same. This as¬ 
sumption will break down if between the epoch of de¬ 
coupling of the inflaton (at time idee) and the onset of 
inflation (at time U) some species annihilates and hence 
causes the temperature of the radiation to be more than 
that of the inflaton at ti. Let us see what are the conse¬ 
quences of relaxing the assumption of Eq. (7). 

We now suppose that at the Planck time, the infla¬ 
ton, BGUT particles, and GUT particles were in thermal 
equilibrium (BGUT refers to Beyond GUT), possibly via 
gravitational interactions. After the Planck scale, at time 
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£dec, the inflaton decouples. A bit later, at a time t\, the 
BGUT particles become non-relativistic and annihilate in 
equilibrium and transfer their entropy to the GUT parti¬ 
cles but not to the inflaton (since it is already decoupled). 
Using b and a for before and after (instantaneous) anni¬ 
hilation, conservation of entropy implies that 


T> f 9*a 'N 3 

Ta \ 9*b ) 


(15) 


where 3*6 is the number of relativistic degrees of freedom 
before the annihilation and g* a is the number of relativis¬ 
tic degrees after the annihilation. 

Assuming no further significant entropy production 
from t\ till ti, the temperature of the inflaton quanta 
at the onset of inflation is 


T i =Tb{a 1 /a i ) , (16) 

while the temperature of the radiation evolves to 

Ty = Taifli/ai) (17) 

1 

= T b (—Y (or/fli) (18) 

\9*a J 

and g le \ = g* a . Then 

Ti = ( T b /Ta )T 7 , (19) 


We thus conclude that any scenario of just enough in¬ 
flation, in which the pre-inflationary universe was radi¬ 
ation dominated and in which the inflaton quanta too 
were in thermal equilibrium at some early epoch, re¬ 
quires an extremely large number of relativistic species 
to be present in the pre-inflationary radiation plasma to 
be consistent with the CMB data. 


IV. DISCUSSION AND CONCLUSION 


We now present a physical interpretation of our re¬ 
sults. The coth[fc/(2T)] factor in our power spectrum 
leads to an enhancement of power on large scales. To 
suppress this requires putting an upper bound on the 
comoving temperature T. If we express the coth fac¬ 
tor for the current horizon scale as coth[fco/(2a*7K)] = 
coth[iJ inf /(2T*)], then the upper bound on T translates 
into an upper bound on T*. For a small value of SN , 
i.e., if the duration of inflation is just enough to solve 
the horizon and flatness problems, Ti ft %■ Now if the 
(frozen) inflaton temperature at the beginning of infla¬ 
tion is the same as that of the thermal radiation in the 
pre-inflationary radiation era, then one gets an upper 
bound on the radiation temperature 


Ty < 4 X 10 13 


GeV 


gjnf 

9 x 10 13 GeV ' 


(24) 


and using Eq. (15) and the fact that p$ = p 7 at the 
beginning of inflation we get 


T = 


{ 9*a 
\ 9*b 



VH in {M P . 


( 20 ) 


Combining Eq. (20) with Eqs. (3) and (4), and using 
the bounds in Eqs. (9) and (12), and setting SN = 0.05 
and Hq ft (4400Mpc) -1 gives 

3*6 > 2 x 10 8 3^ 4 ft 10 9 . (21) 


Such a large change in 3 * can not be obtained from 
the annihilation of a single species. However if multi¬ 
ple species annihilate in equilibrium at and 

transfer their entropy to radiation, then 

T = Tbi(ai/a,i), (22) 



= 7ii(—V fa/ai), (23) 

\ 9*anJ 

as in Eq. (18) with g* an = g leh (g* b m+i = g*am-) Then 
we obtain the same result as in Eq. (21) for 3 * 61 , i.e. 
the number of relativistic degrees of freedom before the 
annihilations must be greater than 10 9 . 


To satisfy this low bound on the pre-inflationary radia¬ 
tion temperature and the condition p r = Pc/> = 3H? ni Mp 
at the beginning of inflation then requires g Te \ to be larger 
than 10 11 . Alternatively, the radiation temperature at 
the beginning of inflation can be much larger than the 
inflaton temperature if entropy was released in the Uni¬ 
verse after the inflaton decoupled. Then while T will be 
less than 4 x 10 13 GeV, T 7 can be as high as 1 x 10 16 GeV 
with 3 re 1 ~ 100. But then 3 * before entropy release is 

(Ty/Tfgrel ~ 10 9 . 

In models of thermal inflation invoked to solve the 
moduli problem [24-26] primordial inflation is followed 
by a later phase of ‘thermal’ inflation which is driven 
by the potential energy of a ‘flaton’ field trapped at the 
origin at the minimum of its thermal effective potential 
[24, 25, 27]. Density perturbations generated during the 
short later phase of thermal inflation influence the spec¬ 
trum only on very small scales while the larger scale per¬ 
turbations which influence the large scale structure and 
cosmic microwave background are largely generated dur¬ 
ing the earlier primordial inflation [28]. Since the flaton 
has thermal interactions one would expect that the flaton 
quanta would have a thermal Bose-Einstein distribution, 
which has so far not been considered. In models of warm 
inflation a thermal bath is generated by dissipation of the 
inflaton field and is maintained in the Universe during 
inflation. If the inflaton fluctuations are also in thermal 
equilibrium with the thermal bath there w T ill again be a 
coth term in the density perturbation spectrum, as seen, 
for example, in Refs. [29-31]. In both thermal and warm 
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inflation, the presence of the coth term would affect den¬ 
sity perturbations with large power on the relevant large 
scales, as in Ref. [4], and can also lead to further con¬ 
straints as discussed above if one also imposes a finite 
duration of inflation. We plan to study these scenarios 
in greater detail in the future. 

In conclusion, we have explored a scenario of inflation 
preceded by a radiation dominated era. We presume that 
the inflaton too could have been in thermal equilibrium 
with the radiation at some very early era and so we con¬ 
sider a frozen thermal state for the inflaton field with 
a Bose-Einstein distribution for the (decoupled) inflaton 
quanta. There is then an upper bound on the comov¬ 
ing relic temperature for the inflaton [4, 8]. Then, in 
these inflation scenarios and with just enough e-foldings 
to solve the horizon and flatness problem, which is mo¬ 
tivated by constraints on large inflaton field excursions 
and are also indicated by the CMB data, we find a lower 
bound on the Hubble parameter during inflation which 
is in severe conflict with the upper bound on the Hubble 
parameter from CMB observations. To resolve this con¬ 
flict one has to allow for an extremely large number of 
relativistic degrees of freedom in the thermal bath in the 
very early Universe, greater than 10 9 or 10 1 11 depending 


on whether or not there has been entropy transfer to ra¬ 
diation after the inflaton decoupled. We are not aware of 
any scenario that may give such large values of g * in the 
early Universe. 
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